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Abstract 

We calculate the one-particle hadronic spectra and correlation functions of pions based on a hy- 
drodynamical model. Parameters in the model are so chosen that the one-particle spectra reproduce 
experimental results of -^s = ISOAGeV Au-f An collisions at RHIC. Based on the numerical solution, 
we discuss the space-time evolution of the fluid. Two-pion correlation functions are also discussed. 
Our numerical solution suggests the formation of the quark-gluon plasma with large volume and low 
net baryon density. 

1 Introduction 

The new state of matter, the quark-gluon plasma (QGP), is expected to be produced in relativistic 
heavy ion collisions. Though the recent experimental data of ^/s — ISOAGeV Au-I-Au collisions 
at RHIC have already been presented f^, the collision process is too complicated to understand 
clearly due to many-body interactions and multiparticle productions. Therefore we need a sim- 
ple phenomenological description as a basis of the further investigation. In this talk, based on a 
hydrodynamical-model calculation, we try to give a possible explanation for the RHIC experimental 
results of the single-particle spectra and correlation function 



2 Model and Particle Spectra 

Assuming that local thermal and chemical equilibrium are established, the dynamical evolution of the 
system can be described by a relativistic hydrodynamics. We numerically solve the hydrodynamical 
equations for the perfect fluid and baryon number conservation with cylindrical symmetry |3| , where 
we take both longitudinal and transverse dynamical expansion into account. We use an equation of 
state (EoS) with a phase transition of first order Q. 

We put the initial time as r = \/t^ — = 1.0 fm. Initial longitudinal velocity is assumed to be 
Bjorken's solution, Vz = z/t. Initial transverse flow velocity is neglected. We parameterize the initial 
energy density distribution and net baryon number distribution as in Fig. ^. 

For both RHIC and SPS, we adjust the parameters in the model so that calculated single-particle 
spectra reproduce the experimental results. The parameter sets for both experiments are summarized 
in Table. 0. 

Once the freeze-out hypersurface is fixed, one can calculate the single-particle spectra via Cooper- 
Frye formula. We take account of not only hadrons directly emitted from the freeze-out hypersurface 
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for RHIC Au+Au 1 30 AGeV for RHIC Au+Au 1 30 AGeV 




Figure 1: Left: Initial energy density distribution. Right: net baryon number distribution. 



Table 1: Parameters and output from the fluid. 





Au+Au at 130 AGeV 


Pb+Pb at 17.4 AGeV 


-»-/max 


6.0 GeV/fm^ 


5.74 GeV/fm^ 


^max 


0.125 fm"3 


0.7 fm-3 


Total riB 


132 


305 


in 


125 MeV 


123 MeV 




76.1 MeV 


280.2 MeV 


("t)|,7|<0.1 


0.51 


0.467 



but also the ones produced in the decay of resonances. Note that the resonances are also emitted 
from the freeze-out hypersurface. We include decay processes p 2n, uj 3n, rj 3n, K* nK, 
and A Nn. 

Figure ^ shows pseudorapidity distribution of charged hadrons and transverse momentum distri- 
bution of negatively charged hadrons. Preliminary experimental data in the pseudorapidity distri- 
bution are taken from the PHOBOS Collaboration |5j . Experimental data in the distribution are 
taken from the STAR Collaboration B. Our results well reproduce the experimental data. 




Pseudorapidity 7]p k-r [GeV/r] 

Figure 2: Single-particle spectra at RHIC. Left: Pseudorapidity distribution. Right: Transverse momen- 
tum distribution. In both figures, solid lines, dashed lines and symbols denote total number of particles 
(i.e., including resonance decay), direct components and experimental data, respectively. 

The parameter set in Table. ^ tells us that the energy density is not extremely high but sufficient 
for the QGP production. The value 6.0 GeV/fm^ is only 5% higher than the one of SPS, 5.74 
GeV/fm^ at the same initial time, r = 1.0 fm. Our value is smaller than different calculations M. 
This is due to the different thermalization time and the different initial energy density profile. Major 
discrepancy between the RHIC and the SPS can be seen in the net baryon density and longitudinal 
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extension of the fluid. The net baryon number at the RHIC is much smaUer than the one at the SPS. 
Longitudinal extension along jy-axis at the RHIC is twice as large as at the SPS. Hence, much higher 
colliding energy at the RHIC is converted into the large volume of the hot matter rather than high 
energy density. We can also see that the transverse flow vt ~ tanhYr becomes stronger at RHIC 
than at SPS. 

It is well known as Hanbury-Brown Twiss effect that one can get the information on the source 
of particles through two-particle correlation function. Based on the numerical solution of the hy- 
drodynamical equation, we calculate the two-pion correlation functions. Here we consider only the 
pions directly emitted from the hypersurface, for simplicity. The HBT radii are extracted from the 
correlation function via a Gaussian fltting function 



C2fit(K, g) = 1 + Aexp [-Rl^,{K)q. 
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Figure 3: HBT radii at RHIC. From left to right, i?sidc,-Rout and i?iong, respectively. Experimental data 
are take from the STAR Collaboration [^. 

Transverse mass dependence of the HBT radii is displayed in Fig. ^. The sideward HBT radii, 
7?aide, show almost consistent result with the experiment. The outward and longitudinal HBT radii 
are larger than the experimental data. The outward HBT radii show different behavior with respect 
to Mt- The experimental result strongly depends on Mt but our result does not show the strong 
dependence on Mt- We also show the i?out/i?sidc which should be a good indicator of the long-lived 
fluid Q in Fig. ^ Our calculation agrees with SPS data, in which the ratio monotonically increases 
with Kt- As shown in Table. |l|. Our numerical solution of the hydrodynamical equations does not 
show major difference in lifetime between the RHIC and the SPS because the initial energy density 
at RHIC is only slightly higher than that at SPS. In fact, the time duration at RHIC is about 6 fm. 
Hence, the ratio shows similar behavior both qualitatively and quantitatively in our result. However, 
the experimental result at RHIC clearly decreases with Mt and the value is smaller than our result. 
For the definite conclusion of this puzzle, we need more and more theoretical investigation. 
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Figure 4: i?out/-Rsido. Left: RHIC. Right: SPS. 
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3 Summary 



In summary, we analyze the one-particle distribution and HBT radii based on a hydrodynamical 
model. Our model successfully reproduce the one-particle distribution at RHIC experiment. The 
numerical solution shows the QGP production with low net baryon number density and large volume. 
The maximum energy density is not much higher than the SPS case. HBT radii are also investigated. 
Our result in Fig. ^shows roughly agreement with the experiment. Figure ^ may indicate different 
particle emission mechanism at RHIC. 
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